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Abstract 

We consider square matrices A that commute with a fixed square matrix K, both 
with entries in a field F not of characteristic 2. When K 2 = I, Tao and Yasuda 
defined A to be generalized centrosymmetric with respect to K. When K 2 = —I, 
we define A to be pseudo-centrosymmetric with respect to K; we show that the 
determinant of every even-order pseudo-centrosymmetric matrix is the sum of two 
squares over F, as long as —1 is not a square in F. When a pseudo-centrosymmetric 
matrix A contains only integral entries and is pseudo-centrosymmetric with respect 
to a matrix with rational entries, the determinant of A is the sum of two inte- 
gral squares. This result, when specialized to when K is the even-order alternating 
exchange matrix, applies to enumerative combinatorics. Using solely matrix-based 
methods, we reprove a weak form of Jockusch's theorem for enumerating perfect 
matchings of 2-even symmetric graphs. As a corollary, we reprove that the number 
of domino tilings of regions known as Aztec diamonds and Aztec pillows is a sum 
of two integral squares. 

Key words: centrosymmetric, anti-involutory, pseudo-centrosymmetric, 
determinant, alternating centrosymmetric matrix, Jockusch, 2-even symmetric 
graph, Kasteleyn-Percus, domino tiling, Aztec diamond, Aztec pillow 



1 Introduction 



This article is divided into two halves; results concerning certain types of 
matrices from the first half apply to enumerative combinatorics in the second 
half. An outline of the structure of the article follows. 

Centrosymmetric matrices have been studied in great detail; in Section [2TTT we 
use their definition as a starting point to define an extension we call pseudo- 
centrosymmetric matrices, defined over a field F, not of characteristic 2. The- 
orem[2]proves that a subclass of these matrices have determinants that are the 
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sum of two squares in F. Specializing further in Section [2721 we show that when 
the matrix is alternating centro symmetric, its determinant has a nice symmet- 
ric form. Returning to the general case in Section 12.31 Theorem El establishes 
that when —1 is not a square in F, even-order pseudo-centrosymmetric ma- 
trices over F have determinants that are a sum of two squares in F. A direct 
corollary is that when the entries are integers, the determinant is a sum of 
integral squares. 

In Section [21 we apply these results on matrices to the question of count- 
ing perfect matchings of graphs. Jockusch proved that the number of perfect 
matchings of a 2-even symmetric graph is always the sum of squares. In Sec- 
tion [37TJ we are able to reprove a weaker version of Jockusch's theorem using 
only matrix methods in the case when the graph in question is embedded in 
the square lattice and the center of rotational symmetry is in the center of a 
unit lattice square. (A restriction of this type is necessary, as shown in Section 
13.31 ) Graphs of this restricted type occur in the study of domino tilings of nice 
regions called Aztec diamonds, as explained in Section [3721 

Future directions of study appear in various remarks throughout the paper. 



2 Matrix-Theoretical Results 

Throughout Sections 12.11 and I2.2[ we define our matrices over an arbitrary 
field F, not of characteristic 2. 

2.1 Pseudo-centrosymmetric matrices with respect to K 

Define J to be the n x n exchange matrix with l's along the cross-diagonal 
(ji^n-i+i) and O's everywhere else. Matrices A such that J A = A J are called 
centro symmetric and have been studied in much detail because of their appli- 
cations in wavelets, partial differential equations, and other areas (see [TH2"] ). 
A matrix such that J A = —AJ is called skew- centro symmetric. Tao and Ya- 
suda define a generalization of these matrices for any involutory matrix K 
[K 2 = I). A matrix A that is centro symmetric with respect to K satisfies 
KA = AK (see [3|TT] ). A matrix A that is skew- centro symmetric with respect 
to K satisfies KA = -AK. 

In the study of generalized Aztec pillows, a related type of matrix arises. 
Define a matrix K to be anti-involutory if K 2 = —I. We will call a ma- 
trix pseudo-centrosymmetric with respect to K if KA = AK, or pseudo-skew- 
centro symmetric with respect to K if KA = —AK. 
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Remark 1 Another definition of centro symmetric matrices is that KAK = 
A. We must be careful with their pseudo-analogues because if KA = AK when 
K 2 = -I, then KAK = -A. 

In this article, we focus on studying even-order pseudo- and pseudo-skew- 
centrosymmetric matrices over F. When K is an n x n matrix for n = 2k 
even, we can write 



K 



Kt K 2 



\K 3 K 4/ 

with each submatrix J*Q being of size k x k. We will explore more general 
anti-involutory matrices K in Section 12.31 but for now we will focus on the 
simple case when Ki = K± = 0. Since K is anti-involutory, K 3 = —K 2 , so 
K is of the form 

K= I . M; 




For such a matrix K, a matrix A that is a pseudo-centrosymmetric with respect 
to K has a simple form for its determinant. 

Theorem 2 Let K and A be matrices defined over F. If K is an anti-involutory 
matrix of the form in Equation (Q, and the 2k x 2k matrix A is pseudo- 
centrosymmetric with respect to K, then A has the form 



A 



B CK 2 
-K 2 l C K 2 7 1 BK 2i 



In addition, det A is a sum of two squares over F. Specifically, if i = 
then over F[i], we have det A = det(-B + iC) det (B — iC). 



PROOF. Calculating the conditions for which a matrix 

/ 

A = 




is pseudo-centrosymmetric with respect to K for matrices K of the form given 
in Equation (TjQ) gives us that K 2 1 A 2 = —A 3 K 2 and A 4 = K 2 1 A 1 K 2 . This 
means we can write 

( B CK 2 

A = 

{-K^C K 2 ~ 1 BK 2/ 

for B = A% and C = A 2 K 2 1 . With this rewriting, det A is simple to compute, 
after adjoining i = y— T to F: 
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det A = det 



B CK 2 
-K 2 X C K 2 l BK 2 



, I 0\ B CK 2 

det | det det 

+iK 2 1 I X-K^C K 2 X BK 2 ) \~iK^ I, 



det 



B-iC CK. 



2 

-It 



\ K 2 \B + iC)K 2/ 
det(B + iC) det(B-iC), 



as desired. 



This is a product of det (5 + iC) and its conjugate, so the determinant of such 
a matrix A is the sum of two squares over F. 

In the case of a matrix A that is pseudo-skew-centrosymmetric with respect 
to K, the analogous result states that A has the form 



.4 



B CK 2 
K 2 X C -K 2 X BK 2 



and that det A = (— l) k det (5 + iC) det (B — iC), which is also a sum of two 
squares over F, possibly up to a sign. The sign term appears because the block 
matrix in the determinant calculation is of the form 

1 B - iC CK 2 



\ 



-K 2 \B + iC)K 2 



with a negative sign appearing from each of the last k rows. Notice that we 
can now explicitly evaluate the determinant of this type of 2k x 2k pseudo- 
and pseudo-skew-centrosymmetric matrix via a smaller k x k determinant. 



In Section [2T3l we will see that the determinants of all even-order pseudo- and 
pseudo-skew-centrosymmetric matrices A can be written as a sum of squares 
(up to a sign) over the base field F, when —1 is not a square in F. 



2.2 Alternating centro symmetric matrices 



We now consider the specific case when the 2k x 2k matrix K is the alter- 
nating exchange matrix — the matrix with its cross- diagonal populated with 
alternating l's and — l's, starting in the upper-right corner. Such a matrix K 
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Fig. 1. The general form of a 6 x 6 alternating centrosymmetric matrix 

is anti-involutory. (Had -ft" been square of odd order, this matrix would have 
been involutory instead of anti-involutory.) 

Definition 3 Let K be the alternating exchange matrix. We define an n x n 
matrix A defined over F to be alternating centrosymmetric with respect to 
K if KA = AK and alternating skew-centrosymmetric with respect to K if 
KA = -AK. 



An equivalent classification of n x n alternating centrosymmetric matrices 
is that their entries satisfy a,ij = (— lY +:j a n +i-i,n+i-j ■ (See Figure [T]). An 
n x n alternating skew-centrosymmetric matrix has entries that satisfy a^j = 

(~^y + ^ +la n+l-i,n+l-j- 

By Theorem [21 we know that the determinant of an alternating centrosym- 
metric matrix is the sum of two squares in F. We now present a different 
version of its determinant with additional symmetry built in. 

We first define a set of fc-member subsets I of [2k] := {1,2,..., 2k}. For any 
subset / of [k], create I by taking I U where i G V if 2k + 1 — i G [k] \ I. 
In this way, each I has k members. We will call subsets of [2k] of this form 
complementary. We define the sets S, S', T, and X" of /c-member subsets of 
[2k]. If I has / elements, 

if / = mod 4, 
if / = 1 mod 4, 
if I = 2 mod 4, or 
if / = 3 mod 4. 



place / into set < 



S 
T 
S' 

r 



Given a 2k x 2k matrix N, we define M(I) to be the k x k submatrix of iV 
with columns restricted to j G /, and rows restricted to the first k rows of N. 

Theorem 4 The formula for the determinant of an alternating centrosym- 
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metric matrix A satisfies 
(let A= ]Tdet(M(/)) det(M(/)) 



./eS 



+ 



£det(M(/))-£det(M(/)) 
JeT ieT' 



PROOF. Let K 2 be the upper-right k x k submatrix of the alternating ex- 
change matrix K. As in the proof of Theorem [2J we consider F[i] and write 
det A = det(B + iC) det(B-iC), where B = A 1 and C = A 2 K 2 X . Calculating 
det(5 + iC) gives x + iy for some x,y G F; we will calculate x and y directly. 



B + iC 



\k+lA 



01,1 + ^i,2fe «i,2 - «ai,2fe-i • • • ai,k + (-1;™ ri ^i,fc+i 
a 2 ,i + ia 2 ,2k a 2 ,2 - ia 2>2k -i ■ ■ ■ a-2,k + (-l) k+1 ia 2 ,k+i 



\dk,l + ^ a fc,2fc a fc,2 — i a k,2k-l ■ ■ ■ a k,k + 



\k+l 



lO'k^+l j 



Define bj to be the column bj 



(aij, a 2j j, . . . , cik,j) T - By linearity of deter- 
minants, det (B + zC) is the sum of 2 k determinants of matrices M with 
dimensions k x k, where in column j we can choose to place either bj or 
i(— l) J+1 &2fc+i-j- Given any determinant of this form, we can convert it into 
a form where the indices of the columns are increasing: I — i\ < ■ • ■ < i r < 
k + 1/2 < z r+1 < • ■ ■ < i k . Note that I is complementary as defined above. 
When we do this and account for changes of sign by interchanging columns, 
the matrix M becomes 



a hil ... a ltir i(-l) k+1 a hir+1 ... i(-l) fc+1 ai jifc ^ 



M 



a 2,h 



\k+l 



a 2,V+i 



*(-1)* + Vh 



ak,i r *(-!) a*,ir+i • • • a k,i k J 



The determinant of this matrix is l) fc+1 ) fc_r times det(M(/)). In partic- 
ular, matrices such that (k — r) =0 mod 2 contribute to x, while matrices 
such that (k — r) = 1 mod 2 contribute to y. In addition, the value of (k — r) 
mod 4 determines the sign of the contribution to the sum. 

This establishes the theorem. 



A similar approach gives an analogous statement for alternating skew-centro- 
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symmetric matrices — we have instead, 



■l) k detA- 



]Tdet(M(/))-£ det(M(/)) 



ieS' 



-\ 2 



£det(M(I))-£det(M(I)) 

JGT JeT' 



Just as in the analogous statement of Theorem El a (— l) k term appears. 



2.3 Pseudo-centrosymmetric matrices in general 



An anti-involutory matrix K behaves like y— T when multiplying a vector. 
When —1 is not a square in F, such a K allows us to create a very nice basis 
of F 2h , similar to the construction of an almost complex structure on IR 2fc . For 
more background on almost complex structures, see Section 5.2 of Goldberg's 
Curvature and Homology [I] . The following lemma will help us prove a general 
result about even-order pseudo- and pseudo-skew-centrosymmetric matrices. 

Lemma 5 If K is a 2k x 2k anti-involutory matrix with entries in F , a field 
not of characteristic 2 in which — 1 is not a square, then F 2k has a basis of 



the form {vi, 



Kvx, . . . , Kv n } for vectors Vj G F 



2k 



PROOF. We will proceed by induction of a set of independent vectors Si = 
{t>i, . . . , vi, Kvi, . . . , Kvi}. Given any nonzero vector v\ G F 2k , we show that 
Si = {vi,Kvi} is linearly independent. Suppose to the contrary that cv 1 = 
Kv\ for some c G F. then c 2 V\ = K(cvi) = K 2 V\ = —v±, so (c 2 + l)t>i = 0, 
which is impossible since \J — 1 ^ F . This proves the base case. 

Now suppose that Si is a set of 21 independent vectors of the form above, 
with I < k. Since 21 < 2k, there exists a vector vi + \ G F 2k not in the span 
of Si. Now we wish to show that Si + i = Si U {vi + i, Kvi + \} is linearly inde- 
pendent. Consider x = v t+1 - proj 5; (i; m ) and y = Kv l+1 - proj 5; the 
components of Vi + i and Kvi + i orthogonal to Si. Notice that Kx = y, so by a 
similar argument to above, x and y are linearly independent. Therefore S^+i 
is linearly independent, and the lemma follows by induction. 

This lemma is the starting-off point for the following theorem; we rely on the 
simpler result from Theorem El 

Theorem 6 Let K be a 2k x 2k anti-involutory matrix with entries in F , a 
field not of characteristic 2 in which —1 is not a square. Let A be a matrix 
(also with entries from F) which satisfies KA = AK . Then the determinant 
of A is a sum of two squares of elements from F . 
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PROOF. Choose k vectors {vj}^ =1 e F 2k such that the set {vi, . . . , v fc , Kvi, . . . , i^v fc } 

forms a basis for F 2fc , as guaranteed by Lemma [5j Define the 2k x 2A; matrix 

V = [vjj • • • Ivfcji^Vxl • • • Since the column vectors form a basis for F 2k , 

the matrix is invertible. Notice that the product KV = [Kvi\ ■ ■ ■ \Kvk\ — Vi| • • • | — v^] = 

-I h 



VK', where K' 







. Therefore, K = VK'V 



Now, if KA = AK, then VK'V~ X A = AVK'V~ X , so K'{V~ l AV) = iy- x AV)K'. 
We see that K' is of the form in Equation (jTJ, so det(A) = det(^ _1 yll / ) is a 
sum of two squares of elements from F, by Theorem [2J 

When A is a pseudo-skew-centrosymmetric matrix, the argument is the same 
up until the last sentence, where we recognize that if KA = —AK, (V~ l AV) 
is pseudo-skew-centrosymmetric with respect to K' , and the remark after The- 
orem [2] implies that the determinant of A is a sum of two squares, possibly up 
to a sign. 

Corollary 7 Let K be an anti-involutory matrix over Q and let A be a matrix 
with integral entries satisfying AK = KA (or AK = —KA). Then det A is a 
sum of two integral squares (possibly up to a sign). 



PROOF. It is clear that det A is an integer. By Theorem O with F = Q, det A 
is the sum of two rational squares. (Possibly up to a sign if AK = —KA.) 
An integer is a sum of two rational squares only if it is a sum of two integral 
squares. This follows from the characterization that an integer n can be written 
as a sum of two squares if and only if its prime factorization contains only 
even powers of primes of the form 4k + 3 (See [3], p. 116). If we can write 
(Pi/Qi) 2 + {Vilq.2? = n for integers p 1 , q u p 2 , and q 2l then p\ + p\ = q\q\n\ 
hence q\q\n is a sum of two integral squares. This implies that q\q\n (and n 
itself) only contain primes of the form 4k + 3 to even powers. 



3 Applications to Enumerative Combinatorics 



The matrix-theoretical results from Sections 12.11 and 12.21 allow us to reprove 
some results from matching theory. 



3.1 Revisiting Jockusch's theorem 



A certain symmetry property of a bipartite graph allows us to say something 
about its number of perfect matchings. 
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Fig. 2. Two graphs that are 2-even symmetric, and two graphs that are not. 

Definition 8 A 2-even symmetric graph G is a connected planar bipartite 
graph such that a 180 degree rotation i? 2 about the origin maps G to itself and 
the length of a path between v and R^iv) is even. 

In particular, if the graph is embedded in a square grid, if you rotate the graph 
by 180 degrees, you get the same graph back, AND the center of rotation is 
not the center of an edge. (See Figure EJ) In [6], William Jockusch proves 
that if a graph is 2-even symmetric then the number of perfect matchings of 
the graph is a sum of squares. Jockusch's result produces a weighted labeling 
function u of the quotient graph G2 involving complex numbers. Once one 
counts the number of weighted matchings associated to u, denoted M U (G 2 ), 
the number of perfect matchings of G is M u (G 2 )M tt (G 2 ), resulting in a sum 
of two squares. 

We can use Theorem [2] to reprove Jockusch's theorem for a subset of 2-even 
symmetric graphs that relies only on the structure of the Kasteleyn-Percus 
matrix of the region, to be discussed shortly hereafter. This subset of graphs 
occurs when the 2-even symmetric graph can be embedded in the square lattice 
with the center of rotation in the center of a square. When we restrict our 
graphs to this type, we can represent the square lattice in the standard x-y 
coordinates by placing vertices at (2k + 1,2/ + 1) for k, I 6 Z, so that the 
center of rotation (0, 0) is the centroid of some square in the lattice. We color 
the vertices white if k + / is even and black if k + / is odd. In Section 13. 3[ we 
present an example that shows that a restriction of this type is necessary. 

Recall that a bipartite graph G = (V, E) is a graph where the vertex set is 
partitioned into two subsets, the "black" and "white" vertices, where there 
are no edges connecting vertices of the same color. Notice that if there is to be 
a perfect matching (or pairing of all vertices using graph edges), there must 
be the same number of black vertices as white vertices. A well-known method 
of counting perfect matchings of a planar bipartite graph relies on taking the 
determinant of a Kasteleyn-Percus matrix. (See [?PD] for background.) The 
definition we will use in this article is that the Kasteleyn-Percus matrix A of a 
bipartite graph G = (V, E) has \ V\/2 rows representing the white vertices and 
\V\/2 columns representing the black vertices. The non-zero entries of A 
are exactly those that have an edge between white vertex Wi and black vertex 
by These entries are all +1 or —1 depending on the position of the edges they 
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Fig. 3. The canonical orientation of edges on the square lattice 

represent in the graph — the restriction is that elementary cycles have a net —1 
product. In the case of the square lattice above, we can satisfy this condition 
easily by giving matrix entries the value —1 if they correspond to edges that 
are of the form e = {v\ v 2) with v\ = (2k — 1,2/ + 1) and v 2 = (2k + 1,2/ + 1) 
and such that v\ is black. This is most easily understood by giving orientations 
to the edges of the lattice as in Figure El and assigning an edge the value +1 
if the edge goes from black to white and the value —1 if the edge goes from 
white to black. The absolute value of the determinant of this matrix counts 
the number of perfect matchings of G. 

With this definition of the Kasteleyn-Percus matrix, we can formulate the 
following theorem. 

Theorem 9 Let G be a 2- even symmetric graph embedded in the square lattice 
with the rotation axis in the center of one square. Under a suitable ordering 
of the vertices, the Kasteleyn-Percus matrix A of G is alternating centrosym- 
metric. 

PROOF. We label the black and white vertices to determine the positions 
of the +1 and —1 entries in A. After an initial labeling, we interchange rows 
and columns as necessary to manipulate the matrix into being alternating 
centrosymmetric, as follows. 

Embedded in this lattice, half the vertices of G lie above the horizontal line 
through the origin. Coloring the vertices of G the color they inherit from the 
lattice coloring above, R2 takes vertices to counterparts of the same color so 
for some m, we have m vertices of each color in the upper half of the graph 
and 4m vertices in all. Label all white vertices v in the upper half of graph 
with values 1 to m, and do the same for black vertices w. For each vertex x 
with value i, label Rz(x) with value 2m + 1 — %. (See FigureHJfor an example.) 
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1 1-1 
1 1 
110 10 

1 1-10 
1-10 

1 1-10 

Fig. 4. Example of the labeling procedure in the proof of Theorem [9j 

From this initial numbering of vertices, we wish to modify some labels so that 
the labels in each row of the square lattice are of the same parity. Note that 
each vertex and its counterpart have opposite parity. Start with the top row. 
For each vertex that is labeled with an even number, switch its label with its 
counterpart. In this way, all elements in the top row will have an odd value. 
For the second row, exchange a vertex's label with its counterpart's label if it 
has an odd value. Continue in this fashion until all odd rows have vertices with 
odd labels and all even rows have vertices with even values. After determining 
all rows above the horizontal line through the origin, the rest of the rows come 
for free. 

By this construction, for any horizontal edge (vi,Wj) in G', we know that % 
and j are of the same parity. Similarly, we know that any vertical edge (vi, Wj) 
has % and j of opposite parity. In addition, this implies that the rotation of a 
vertical edge by R 2 results in the opposite sign appearing in the Kasteleyn- 
Percus matrix. 

Since (v^Wj) is a horizontal edge if and only if (t>2 m+ i_j, u^m+i-j) is a hor- 
izontal edge, a +1 appears in position an^ for i + j even if and only if 
a(2m+i-i,2m+i-j) is +1. Similarly, if (v h Wj) is a vertical edge, then so is (v 2m +i-i, w 2m +i-j), 
and their entries in A are opposite. This occurs exactly where i+j is odd. All 
other entries are zero, so for those a^j we have a^j) = ±a(2 m +i-i,2m+i-j) ■ 

These conditions imply that entry a^j) equals (— l)* + - ? a(2 m +i-i ) 2m+i-j), which 
implies A is alternating centrosymmetric, as desired. 

With this theorem and Corollary [71 we have the following corollary. 

Corollary 10 The number of perfect matchings of a 2- even- symmetric graph 
embedded in the square lattice with the center of rotation in the center of one 
unit square is a sum of two integral squares. 

Remark 11 This theorem is weaker than Jockusch's original theorem as it 
applies to fewer regions. However, it allow us to prove a sum-of- squares result 
using only matrix-based methods. 
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Fig. 5. Examples of an Aztec diamond, an Aztec pillow, and a generalized Aztec 
pillow 




Fig. 6. The dual graphs to the regions in Figure [5j 
3.2 Applications to generalized Aztec pillows 

An application of Corollary [10] has to do with domino tilings of rotationally- 
symmetric regions made up of unit squares. A domino tiling of a region is a 
complete covering of the region with non-overlapping 2x1 and 1x2 rectangles 
(or dominoes). We can associate to any region its dual graph, a graph with 
vertices representing the unit squares and edges between vertices representing 
adjacent squares. In effect, counting the number of domino tilings of a given 
region is the same as counting the number of perfect matchings of the corre- 
sponding dual graph. Many results in this vein are presented in James Propp's 
survey article [8]. Some regions fit the framework from the previous section 
particularly well; examples are pictured in Figure [5] and their corresponding 
dual graphs are pictured in Figure [6j 

An Aztec diamond is the union of the 2n(n + 1) unit squares with integral 
vertices (x, y) such that |x| + \y\ < n + 1. Aztec pillows were introduced in 
[SJ and explored more in depth in [SJ. They are also rotationally-symmetric 
regions composed of unit squares with their "steps" along the northwest and 
southeast diagonals having height one and length three. Lastly, we include 
generalized Aztec pillows, where all steps off the central band of squares are 
of height one and odd length. Notice that this implies that Aztec diamonds 
and regular Aztec pillows are also generalized Aztec pillows. While generalized 
Aztec pillows need not be rotationally symmetric, when we restrict to those 
that are rotationally symmetric, we now have many regions whose dual graphs 
are 2-even-symmetric graphs, so we have the following corollary. 

Corollary 12 The number of domino tilings of any rotationally-symmetric 
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Fig. 7. A 2-even symmetric graph whose Kasteleyn-Percus matrix is not alternating 
centr osymmetr ic . 

generalized Aztec pillow, which includes all Aztec diamonds and Aztec pillows, 
is a sum of two integral squares. 



3.3 Problems extending Theorem^ 



Unfortunately, the methods from Section 13.11 do not allow us to completely 
reprove Jockusch's theorem, as we highlight with the following example. Con- 
sider the graph in Figure This graph is clearly bipartite, can be embedded 
in the square grid, and is 2-even symmetric. However, the following theorem 
holds. 

Theorem 13 No reordering of the vertices transforms the Kasteleyn-Percus 
matrix of the graph in Figure [?| into an alternating centro symmetric matrix. 



PROOF. Notice that our relabeling trick from Section I3TT1 will not work here, 
as the would contain both even- and odd-labeled vertices. In addition, 

any rearrangement of the vertices could produce neither an alternating cen- 
trosymmetric matrix nor an alternating skew-centrosymmetric matrix. We can 
see this as follows. Vertex a is adjacent to vertices 1 and 2, while vertex c is 
adjacent to vertices 3 and 4, all by +1- weighted edges through our scheme 
from Figure [3J This implies that whatever rows of our rearranged matrix A 
the vertices a and c are in, the entries used in the rows will be complementary 
(such as of type (I) 1010 and 0101, type (II) 1001 and 0110, or type (III) 1100 
and 0011). This is true also with b and d, with the addition of some signs. 
Since vertices b and d share one vertex each with a and c, they can not have 
complementary pairs of the same type (I, II, or III). In order for the non-zero 
entries of A to match up correctly, a's and c's rows must be the center two rows 
or the first and last rows. Therefore only types (I) and (III) are valid types, 
and there must be one of each for the two pairs. Unfortunately, this can not 
possibly work when considering the sign conventions necessary for a matrix 
to be alternating centrosymmetric or alternating skew-centrosymmetric. 
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Remark 14 As alluded to in the introduction of Section 13. 1\ there are multi- 
ple definitions of a Kasteleyn-Percus matrix. Perhaps it is possible to reprove 
Jockusch's theorem using another matrix interpretation. 

Remark 15 Although we may not be able to reprove Jockusch's theorem using 
matrix methods for all 2- even symmetric graphs, perhaps the condition that the 
graph be embedded in the square lattice can be relaxed. 
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